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ABSTRACT 

In this paper we have applied a new approximation technique to solve non-linear partial differential equation. 
The approximation technique is called Homotopy perturbation method (HPM). The main difference between traditional 
perturbation method and this one is that it can be applied for even higher values of the parameter, where as traditional 
perturbation method can be applied only for lower values of parameter. It means that when the value of the parameter is 
less than one only then the traditional perturbation can be applied. In this paper we have considered a highly non-linear 
partial differential equation and found the approximate solution using HPM for two types of initial conditions. Then we 
have drawn two-dimensional and three dimensional graphs from the solution of the equations which demonstrates the 
physical situation of the solution for different values of the parameter. This gives us the clear picture of the range of the 
variables for which the normal solution exists and for what values of the variable the chaotic situation arises. 
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1. INTRODUCTION 

Chaos or chaotic system received great attention among mathematicians and physicists. Because it stemps out 
from natural phenomena. Mathematically one can get this studying linear or non-linear dynamical system or non-linear 
partial differential equations. Unfortunately to get a close form of the solution of any of them is almost beyond our present 
knowledge. However, one can use some kind of analytical or numerical methods to get an approximate solution. 

Recently, J. H. He [1, 2] found an ingenious method to solve non-linear differential equation, which is called 
"Homotopy Perturbation Method" (HPM). He applied this method to solve different types of non-linear equations, such as 
Van der Pol- duffing equations [3] for different types of oscillator [4, 5, 6] problems and so on. 

Here we have applied HPM technique to find the analytical solution of a non-linear differential equation with 
different types of initial conditions. Taking different values of the parameters and for a wide range of time we have found 
the graphical representation of the analytical solution for all domain of response. Following He's paper [2] we have 
presented the basic prescription of the method in section 2. In section 3 we have applied this technique to solve non-linear 
partial differential equation with two different types of initial conditions and found an approximate solution of the problem. 
In section 4 we have given the graphical representation of the solution for two different initial conditions. In last section we 
have discussed the nature of the solution for different values of the parameters and range of variables 

Homotopy Perturbation Prescription 

J. H. He [1] gave a prescription to solve nonlinear differential equation which can be expressed in the following 

way. 
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Let us consider a nonlinear differential equation of the form 

A(w)-/(r) = 0, reQ (l) 

Subject to the boundary conditions: 

B(u,du/dn) = 0, reT (2) 

where A is a general differential operator, B a boundary operator, f(l") a known analytical function and T is 
the boundary of the domain Q . 

In general one can divide the operator A into two parts: linear and non- linear. That means 

A = L + N 

where L is linear and N is non-linear. 
Hence, equation (1) can now be rewritten as 

L(u) + N(u) -f(r) = 0, r e Q (3) 

By the homotopy technique, one can construct a homotopy in the following way v(r, p) :fix [0,1] ~~ ► R 
which satisfies 

H(v, p) = (l- p)[L(v) - L(u )] + p[A(v) - /(r)] = 0, p e [0,1] , r e Q (4) 
or (v, p) = L(v) - L( Mo ) + pL(u ) + p[N(v) - f(r)] = (5) 

Where p G [0,1] is an embedding parameter, U Q is an initial approximation of equation (1) which satisfies the 
boundary conditions. 

Obviously, from equations (4) and (5) we will have: 

H(v,0) = L(v)-L(u ) = (6) 
tf(v,l) = A(v)-/(r) = (7) 

The changing process of p from zero to unity is just that of v(r, p) from U Q (r) to ll(r) . In topology, this is 

called deformation, and L(v) — L(u ) and A(v) — f(f) are called homotopy. According to the HPM, we can first use 
the embedding parameter p as a "small parameter", and assume that the solution of equations (4) and (5) can be written as 
a power series in p 

v = v Q +pv [ +p 2 v 2 + (8) 

Setting p = 1 results in the approximate solution of equation (1): 

U = lim V = V + V, + V 2 + (9) 
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The combination of the perturbation method and the homotopy method is called the homotopy perturbation 
method (HPM), which has eliminated the limitations of the traditional perturbation methods. On the other hand, this 
technique can have full advantage of the traditional perturbation techniques. The series (9) is convergent for most cases. 

2. SOLUTION OF NON-LINEAR PARTIAL DIFFERENTIAL EQUATION 

In this section we will show how He's homotopy method can be applied non-linear partial differential equation 
with different types of initial conditions. In this view we have considered a highly non-linear partial differential equation 
with two types of initial conditions, and used He's homotopy method to solve these problems. The two different types of 
initial conditions considered here is just to see the sensitivity of the behaviour of the solutions for different types of initial 
conditions. Because in any dynamical system it is very important to check whether this types of equations is sensitive to 
initial conditions. 

3.1. Solution of Non-Linear Equation with Initial Conditions 

Let us consider a highly non-linear partial differential equation 

U tt +aU xx +pU + yU'=0 (10) 

and try to find the solution with the following initial conditions: 

U(X,0) = X 3 and U,(X,0)=0 (11) 
In order to solve equation (10) using HPM, a homotopy- perturbation method can be constructed as follows 

d 2 v d 2 v d 2 v 

H(v,p) = (l-p) — + p(— + a — + j3v + yv 3 ) = (12) 
ot ot ox 

Substituting U = Jim V = lim[v + p V, + p 2 V, + into equation (12) and rearranging the resultant 

p->i p->i 

equation based on powers of p -terms, we can find the following equations: 



i d 2 v, d 2 v,, 

P :-^- + a—f + fiv +yv = (14) 
ot ox 



2 d 2 v 2 d 2 v 1 2 n 

p :— -r + a — l - + j3v l +yv v, =0 (15) 
ot ox 

Solving equation (13) we get 

— - = A (16) 
dt 

v =At + B (17) 
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Before using the initial condition in (16) and (17) let us clarify one important point. If we back to earlier 
substitution that is 

u(x, t) = lim v = v + Vj + v 2 + 

P ->i 

and taking the initial conditions : U (x,0) = X 3 and U t (x,0) = 

which implies, V = X 3 and V l = V 2 = = 

Using the initial conditions (1 1), at t = we get A = and B = X 3 
Hence, finally we get the solution 

V Q =X 3 (18) 
From this equation we find 

-^ = 3x 2 , — f = 6x 
dx dt 2 



dv d v 

After substituting the value of V , and — in (14) 

dx dt 

we get the equation for V l 

d 2 v, _ n 3 9 , 



L = -(6«x + ^x ) 



a^ 2 

Integrating 

— L = -(6 ax + fix +yx )t + C 
dt 

Again integrating we arrived at 

v 1 =-^(6ax + J3x 3 +yx 9 )t 2 +Ct + D 

Using the initial conditions in the above equations we get 

at t = 0, ^ = = C, v 1 (x,0) = = D 
dt 

Putting the values of C and D we get 

v l =-^(6ax + ft x 3 +y x 9 )t 2 (19) 
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which is the solution for Vj 

Now to solve equation (15) for V 2 let us find 

d ^ = --(6a + 3/Jx 2 +9yx s )t 2 
dx 2 



d\ __ 1 



7 \ .2 



(6j3x + 12yx')t 



dx 2 2 

Using these in equation (15) we get 

<3 2 v 7 , d 2 v, _ 2 N 

= Ul2a x + fi 1 x 3 +90a r x 7 +APyx' +3y 2 x")t 1 

Integrating, 

^5 1 

-^ = -(12a/?x + /? 2 x 3 +90a^x 7 +4/?^x 9 + 3^V 5 )r 3 +C 
3f 6 

Again integrating, 

v 2 =^-(12a/?x + /? 2 x 3 +90a^x 7 + 4j3yx 9 +3y 2 x 15 )t 4 +C't + D' 

Again the initial conditions (11) implies that 

C'=0 and D'=0 
Therefore 

v 2 =^-(12a/?x + /? 2 x 3 +90a^x 7 +4/?^x 9 +3^ 2 x 15 )r 4 (20) 

According to HPM we can write the solution of (10) up to second order of p as : 

u(x,t) = lim v = lim[v + pv { + p 2 v 2 ] 

Setting p = 1 the above equation becomes, 

u(x,t) = V + Vj + v 2 

Substituting the values of V , V { and V 2 

u(x,t) = — [24x 3 -(72ax + l2j3x 3 + 12yx 9 )t 2 
24 
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+ (l2aj3x + j3 2 x 3 +90ayx 7 +4j3yx 9 +3y 2 x 15 )t 4 ] (21) 

This is the approximate solution of the non-linear partial differential equation (10) with initial conditions (11). 
3.2. Solution of Same Non-Linear Differential Equation with another Initial Conditions 

Again considering the same differential equation 

U tt +aU xx + j3U + yU 3 =0 (22) 

with different type of initial conditions: 

U(X,0) = and U t O,0) = e Xx (23) 
Following the same procedure for this problem we get, 

U(x,t) = ^e Ax [\20t-(20aZ 2 + 20 J3)t 3 + (a 2 Z 4 +2a/U 2 + J3 2 )t 5 ] 

— — e 3Ax [42yt 5 -(\9ayl 2 +UBy)t 7 ] + — y 2 e 5Ax t 9 (24) 
840 480 

This is the approximate analytical solution of the partial differential equation (22) with another type of initial 
conditions. 

With the help of HPM we are able to solve highly nonlinear higher order partial differential equation with 
different types of initial conditions. 

In the following sections we will find the graphical representation of the solutions of non-linear differential 
equations (10) and (22) found in (21) and (24) for different values of parameters (X, |3, y and for a wide range of x and 
t values. From which we can draw some conclusion of the results. 

4. GRAPHICAL REPRESENTATION OF THE SOLUTIONS OF (10) AND (22) 

In this section we have drawn three-dimensional and two-dimensional graphs from the solution (21) of non-linear 
differential equation (10) with the given initial conditions taking larger and smaller values of the parameter and for a wide 
range of the values of the variables x and t . After a close look to the graphs we have divided the shape of the graphs in 
four kinds. The graphs of similar kind for different parametric values but for a particular range of x and t are depicted 
below. 
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Figure 1 Figure 1' 

Figure 1: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , P = 1 and 
y = l and Values of the Variables x e (0,2.5) and t e (0,.2) . Figure 1' Shows the Corresponding 

Two-Dimensional Graph 
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Figure 2 Figure 2' 

Figure 2: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = . 5 , P = . 6 
and y = .7 and Values of the Variables x e (0, 2 . 5) and t e (0, . 2) . Figure 2' Shows the 
Corresponding Two-Dimensional Graph 




70 












60 
50 












40 
30 
20 
10 














0.5 


1 


1.5 


2 


2.5 



Figure 3 Figure 3' 

Figure 3: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , P = 2 and 
y = 3 and Values of the Variables x e (0,2.5) and t e (0,.2) . Figure 3' Shows the Corresponding 

Two-Dimensional Graph 
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Figure 4 

Figure 4: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 5 , P = 1 and 
y = 2 and Values of the Variables x e (0,2.5) and t e (0,.2) . Figure 4' Shows the Corresponding 

Two-Dimensional Graph 




Figure 5 Figure 5' 

Figure 5: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , P = 5 and 

y = 1 and Values of the Variables x e (0,2.5) and t e (0,.2). Figure 5' Shows the Corresponding 

Two-Dimensional Graph 
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Figure 6 Figure 6' 

Figure 6: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , 

P = 1 and y = 1 and Values of the Variables x e (0, 1000) and t e (0, 1000) Figure 6' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 7 Figure 7' 

Figure 7: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = .5 , 
P = .6 and y = .7 and Values of the Variables x e (0, 1000) and t e (0, 1000). Figure T Shows the 

Corresponding Two-Dimensional Graph 
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Figure 8 Figure 8' 

Figure 8: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , P = 2 and 

y = 3 and Values of the Variables x e (0, 1000) and t e (0, 1000). Figure 8' Shows the 
Corresponding Two-Dimensional Graph 
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Figure 9 Figure 9' 

Figure 9: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , P = 5 and 

y = 1 and Values of the Variables x e (0, 1000) and t e (0, 1000). Figure 9' Shows the 
Corresponding Two-Dimensional Graph 
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Third Kind 
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Figure 10 Figure 10' 

Figure 10: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , p = 1 and 

y = 1 and Values of the Variables x e (700, 1000) and t e (0, 1000). Figure 10' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 11 Figure 11' 

Figure 11: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = .5 , P = .6 and 

y = .7 and Values of the Variables x e (700, 1000) and t e (0, 1000). Figure 11' Shows the 

Corresponding Two-Dimensional Graph 




Figure 12 Figure 12' 

Figure 12: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , P = 2 and 

y = 3 and Values of the Variables x e (700, 1000) and t e (0, 1000). Figure 12' Shows the 

Corresponding Two-Dimensional Graph 
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Fourth Kind 




Figure 13 Figure 13' 

Figure 13: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , 
P = 1 and y = 1 and Values of the Variables x e ( 1 , 1.00001) and t e (0, 1000). 
Figure 13' Shows the Corresponding Two-Dimensional Graph 




Figure 14 Figure 14' 

Figure 14: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = .5 , 
P = .6 and y = -7 and Values of the Variables x e (1, 1.00001) and t e (0, 1000). 
Figure 14' Shows the Corresponding Two-Dimensional Graph 




Figure 15 Figure 15' 

Figure 15: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , 
P = 2 and y = 3 and Values of the Variables x e (1,1.00001) and t e (0, 1000). Figure 15' Shows 

the Corresponding Two-Dimensional Graph 
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Figure 16 Figure 16' 

Figure 16: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = — 1 , 
P = -1 and y = -5 and Values of the Variables x e ( 1 , 1.00001) and t e (0, 1000). Figure 16' 
Shows the Corresponding Two-Dimensional Graph 

Analysis of the Solution of Equation (22) 

The solution of equation (22) is given in equation (24). From this solution we have drawn many three-dimensional 
and two dimensional graphs for different values of the parameters and for a wide range of the variables. The significant 
graphs are shown below. If we look very carefully these graphs, we can see that the shape of the graphs changes for 
different range of variables. A slight change of the graphs has also been seen for positive and negative values of the 
parameters. 

First Kind 




Figure 17 Figure 17' 

Figure 17: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , [3 = 1 , 

y = 1 and A. = 1 and Values of the Variables x e (0, 100) and t e (0, 100). Figure 17' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 18 Figure 18' 

Figure 18: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = — 1 , p = 1 , 

y = 1 and X = 1 and Values of the Variables x e (0, 100) and t e (0, 100). Figure 18' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 19: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = .5 , P = .6 , 

y = .7 and X = .8 and Values of the Variables x e (0, 100) and t e (0, 100). Figure 19' Shows the 

Corresponding Two-Dimensional Graph 

Second Kind 
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Figure 20 Figure 20' 

Figure 20: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , P = 1 , 

y = 1 and X = 1 and Values of the Variables x e (0, 1000) and t e (0, 1000). Figure 20' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 21 Figure 21' 

Figure 21: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = — 1 , 
P = 1 , y = 1 and X = 1 and Values of the Variables x e (0, 1000) and t e (0, 1000). Figure 21' Shows 

the Corresponding Two-Dimensional Graph 

Third Kind 
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Figure 22 Figure 22' 

Figure 22: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = .5 , P = .6 , 

y = .7and A, = .8 and Values of the Variables x e (0,250) and t e (0, 1000). Figure 22' Shows the 

Corresponding Two-Dimensional Graph 

Fourth Kind 
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Figure 23: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = .5 , p = .6 , 

y = .7 and X = 2 and Values of the Variables x e (0, 136) and t e (0, 1000). Figure 23' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 24 

Figure 24: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , P = 1 , 
y = 1 and X = 5 and Values of the Variables x e (0,500) and t e (0, 1000). Figure 24' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 25 Figure 25' 

Figure 25: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , [3 = 1 , 

Y = land X = -1 and Values of the Variables x e (0, 1000) and t e (0, 1000). Figure 25' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 26 Figure 26' 

Figure 26: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , [3 = 1 , 

Y = 1 and X = 1 and Values of the Variables x e (130, 145) and t e (0, 1000). Figure 26' Shows the 

Corresponding Two-Dimensional Graph 
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Seventh Kind 
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Figure 27 Figure 27' 

Figure 27: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 ,(3 = 1 , 

y = 1 and A = 1 and Values of the Variables x e (130, 135) and t e (0, 1000). Figure 27' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 28: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = — 1 , p = 1 , 

y = 1 and A, = 1 and Values of the Variables x e ( 1 30, 145) and t e (0, 1000). Figure 28' Shows the 

Corresponding Two-Dimensional Graph 
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Figure 29: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = — 1 , 
P = 1 , y = 1 and X = 1 and Values of the Variables x e ( 1 30, 1 35) and t e (0, 1 000). Figure 29' 
Shows the Corresponding Two-Dimensional Graph 
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Figure 30 Figure 30' 

Figure 30: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , [3 = 1 , 

y = 1 and 1 = -1 and Values of the Variables x e (130, 140) and t e (0, 1000). Figure 30' Shows 

the Corresponding Two-Dimensional Graph 
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Figure 31: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , P = 1 , 

y = 1 and 1 = -1 and Values of the Variables x e (-135,-12Q and t e (0, 1000). Figure 31' 
Shows the Corresponding Two-Dimensional Graph 

Tenth Kind 




100 110 120 130 140 

Figure 32 Figure 32' 

Figure 32: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , [3 = 1 , 

y = 1 and 1 = 1 and Values of the Variables x e (97, 1000) and t e (0, 1000). Figure 32' Shows 

the Corresponding Two-Dimensional Graph 
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Figure 33: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = .5 , 
P = .6 , y = .7 and X = .8 and Values of the Variables x e (133, 1000) and t e (0, 1000). 
Figure33' Shows the Corresponding Two-Dimensional Graph 
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Figure 34: The Surface Shows the Solution U (x, t) when the Values of the Parameters 
a = .5 , P = .6 , y = .7 and X = 2 and Values of the Variables x e (27, 1000) and t e (0, 1000). 
Figure 34' Shows the Corresponding Two-Dimensional Graph 



Eleventh Kind 
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Figure 35: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = 1 , p = 1 , 

y = 1 and 1 = 1 and Values of the Variables x e (99, 1000) and t e (0, 1000). Figure 35' Shows 

the Corresponding Two-Dimensional Graph 
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Figure 36 Figure 36' 

Figure 36: The Surface Shows the Solution U (x, t) when the Values of the Parameters a = — 1 , 
P = 1 , y = 1 and 1 = 1 and Values of the Variables x e (99, 1000) and t e (0, 1000). Figure 36' 
Shows the Corresponding Two-Dimensional Graph 
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Figure 37 Figure 37' 

Figure 37: The Surface Shows the Solution U (x, t) when the Values of the Parameters 
a = .5 , P = .6 , y = .7 and X = 2 and Values of the Variables x e (28, 1000)and t e (0, 1000). 
Figure 37' Shows the Corresponding Two-Dimensional Graph 

RESULTS AND DISCUSSIONS 

To find the nature of the graphs and to give the physical meaning of the solution of non-linear differential 
equation (10) we have drawn 16 numbers of three dimensional and two dimensional graphs for a wide or close range of 

values of the variables x and t and for different values of the parameters (X, |3, J . These graphs have shown that the 

approximate solution (21) is valid for all values of x and t and for lower and higher parametric values of 01, (3 and J . 

These graphs also indicate that the shape of the graph depends on the values of the variables but not on parametric values. 
Here we can see four different shapes of the graph which are mostly depends on the range of the values of x and t but not 
on any parameter. We can classify these four kinds in the following way. Say for all parametric values but for close range 
of x and t that is if x is from to 2.5 and t is from to 0.2 then let us represent the graph by first kind, which are depicted 
in figures 1 to 5. Similarly, for a wide range values of x and t that means say x is from to 1000 and t is also from 
to 1000 the graphs are same which we call second kind which are depicted in figures 6 to 9. Another shape of the graph 
one can find for close range of x but wide range of t, which we call as third kind. That is the x values say is from 
700 to 1000 but t is from to 1000, which are depicted in figs. 10 to 12. But when the range of x is from 1 to 1.00001 and 
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the range of t is from to 1 or to 1000 we have got another shape of the graphs and we call these types of the shape as 
fourth kind, which are depicted in figures of 13 to 16. The two-dimensional graphs that is for a fixed values of t but for a 
wide range of x the corresponding figures of first kind are depicted in figures from 1' to 5' and that of second kind are 
depicted in figures from 6' to 9' and the third kind are depicted in figures from 10' to 12' and the fourth kind are depicted in 
figures from 13' to 16'. All these graphs shows that the approximate solution that found by HPM method are most accurate 
and there is no chaotic solution in this type of highly non-linear differential equation with this type of initial conditions. 
Also it is clear from the graphs that the solution is valid for lower and higher values of the parameter, which is the main 
advantage of the HPM method. 

The approximate solution of (22) that is for the same non-linear equation with the other type of initial conditions, 
the solution found by HPM method is given in equation (24). To find the nature of the solution we have drawn many three 
dimensional graphs for a wide or short range of values of the variables x and t and for different values of the parameters 
and corresponding two- dimensional graphs for fixed values of t. From these graphs we have chosen only 21 numbers of 
both three-dimensional and two-dimensional graphs to show the differences of shape of the graphs. 

In the first case that is for the differential equation (10) we got only four types of different graphs, which are 
almost same but a minor difference in bend point, which is negligible. But in the second case that is for other type of initial 
conditions we got various types of graphs. This is mainly because of the different type of initial conditions. The main point 
here is that the solution of equation (22) exists for all values of the parameters except negative values of the parameter X 
for a particular range of values of x and t which are depicted in figures from 17 to 24. After certain range of x and t values 
the solution does not exists. And these can be seen for the different range of the values of the x variables and the different 
values of the parameters. This can be easily understood if we look very carefully to the figures: 17 to 24, we can see that 
the solution exists only if x lies between to 136 approximately but for all t and after that there is no solution. This one can 
see from the graphs in figure 35 to figure 37. For negative exponential initial conditions the solution exists in any range of 
values of x and t, which can be seen in figure 25. The solution (24) gives us some irregular shape of the graphs and also 
chaotic solution arises with these initial conditions, whereas in the solution (21) there was no such type of graphs found. 
This is only because of the initial conditions. The chaotic solution found for a very close range of x values but almost for 
all t and even for any parametric value. This one can find in the graphs of figures 26 to 29. Again if we compare 
figure 25, 30, 31 we can see that for - A that is for negative exponential the graph depends on X -values but not on t . 

It is clear from these calculations that the solution is sensitive to the initial conditions. For the same equation due 
to the first type of initial conditions we got regular solution for all x and t. But in the second case we got both regular and 
irregular or chaotic solution. This proves the sensitivity of the initial conditions. From these solutions it is clear that HPM 
is a good semi-analytic approximation method to solve non-linear differential equation or dynamical system. 
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